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We investigate the evolution of cosmological perturbations 
in scenarios with a quintessence scalar field, both analytically 
and numerically. In the tracking regime for quintessence, we 
find the long wavelength solutions for the perturbations of the 
quintessence field. We discuss the possibility of isocurvature 
modes generated by the quintessence sector and their impact 
on observations. 

PACS numbers: 98.80Cq 



I. INTRODUCTION 

Observational data seems to indicate that we live in an 
accelerating Universe Q. As an alternative to the sce- 
nario where the acceleration is fueled by a cosmological 
constant A, models with a scalar field capable of dom- 
inating recently and of developing a negative pressure 
have been proposed . This scalar field which pervades 
our universe has been dubbed quintessence. 

Compared to the cosmological constant, quintessence 
has two important differences. First, quintessence can 
be interpreted as a fluid with a time dependent equation 
of state. Therefore quintessence models may alleviate 
the so-called coincidence problem, which is the appar- 
ent cosmic collusion that the dark energy component is 
fine-tuned in a way that it is starting to dominate the 
energy density of the universe just at the present time. 
And second, in contrast to the cosmological constant, the 
quintessence field can fluctuate J|-Q]. 

One interesting possibility is that the quintessence field 
Q(x,t) can, in combination with the other cosmic fluids 
(radiation, baryons, cold dark matter, etc.), lead not only 
to adiabatic (curvature) perturbations, but to a mixture 
which includes an isocurvature component. Isocurvature 
(or entropy) perturbations appear when the relative en- 
ergy density and pressure perturbations of the different 
fluid species combine to leave the overall curvature per- 
turbations unchanged. In quintessence models, the pres- 
ence of potentially relevant isocurvature modes could be 
generic, just as in multi-field inflationary models Q. In- 
deed, quintessence is constructed is such a way that it 
is an unthermalized component subdominant for most of 
the history of the Universe. Since quintessence is un- 
coupled from the rest of matter because of astrophysi- 
cal and cosmological constraints M , its fluctuations may 



lead to an isocurvature component, whose nature will be 
preserved except for the known integrated feeding of the 
adiabatic component, expressed by the relation po| , |lT|| : 
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where C is the gauge-invariant curvature perturbation, 
a dot denotes time derivative, w = p/p and c 2 = p/p 
denote respectively the equation of state and the speed 
of sound of the total matter content, p na d is the total 
non-adiabatic pressure perturbation and we employ units 
such that 8irG = 1. 

This latter effect was studied [|l2| in the context of 
axion perturbations, and [p3| in the context of the baryon 
isocurvature model [H. This effect is also responsible 
for the growth of super-Hubble adiabatic perturbations 
during preheati ng Rig] . Considering axions as Cold Dark 
Matter (CDM) ]12|, an isocurvature perturbation due to 
the angle misalignment produced during inflation induces 
an adiabatic component of comparable amplitude at the 
moment of reentry of the perturbation inside the Hubble 
radius. This is due to the fact that the CDM component 
is going to dominate about the time of decoupling, and 
thus the integrated effect is almost completed by the time 
that mode reenters inside the Hubble radius. 

However, the quintessence case is different from the ax- 
ion/CDM since in most of the models quintessence fluc- 
tuations are damped inside the Hubble radius. This is 
required in order to minimize the impact of an additional 
dynamical degree of freedom on structure formation. The 
quintessence and the axion/CDM differ also in either one 
of two ways: i) Q was still a negligible component be- 
fore the time of decoupling, or ii) Q was comparable to 
normal matter, Qq = 0(1), but in a so-called tracking 
regime p||l6| whereby its equation of state wq = Pq/pq 
was approximately that of dust or radiation — whichever 
was dominating at the time. In the first case (which hap- 
pens, for example, in the PNGB scenario pJ§ ) isocurva- 
ture perturbations are irrelevant simply because the field 
Q(x,t) is a negligible component until a redshift of at 
least z ~ 10. In the second case the energy density in 
Q need not be small, however due to the tracking of the 
quintessence field, perturbations in the Q-fluid behaved 
similarly to the perturbations in the background for most 
of the observable history of the Universe, and isocurva- 
ture perturbations are therefore suppressed until the end 
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of the tracking phase. In either primordial isocur- 

vature perturbation could still be present, but it would 
not have had enough time to induce an adiabatic com- 
ponent. 

From an observational point of view, isocurvature per- 
turbations have a very distinctive imprint on the spec- 
trum of the temperature anisotropy of the Cosmic Mi- 
crowave Background (CMB) With the accu " 
racy that future CMB experiments such as MAP 
and Planck |2^| will be able to reach, the constraints on 
the ratio of uncorrelated isocurvature perturbations in 
CDM-radiation to the adiabatic component will be of the 
order of percents 23 , 24|] . Recently the impact of generic 
isocurvature modes on the estimation of cosmological pa- 
rameters has been also investigated Therefore it is 
important to understand the evolution of isocurvature 
modes in quintessence models where an unthermalized, 
uncoupled relic survives until the present era, and domi- 
nates very recently. In most of the literature a primordial 
adiabatic spectrum for quintessence perturbations is as- 
sumed. If the notion of adiabaticity among different com- 
ponents is related to their thermal equilibrium, then the 
weakly coupled nature of quintessence could evade this 
condition. The primordial spectrum could be generated 
during inflation and/or influenced through its evolution 
until the decoupling time. 

The outline of the paper is as follows. In section II 
we give necessary notions of the background evolution 
of quintessence models. In section III we study the evo- 
lution of cosmological perturbations and we identify the 
attractor solution for quintessence perturbations during 
the tracking regime. In section IV we study the evolution 
of isocurvature perturbations and their feedback on the 
adiabatic component. In section V we discuss the initial 
condition for quintessence perturbations after nucleosyn- 
thesys. We conclude in section VI. 



II. BACKGROUND EVOLUTION WITH 
QUINTESSENCE 

For simplicity we consider only radiation, pressureless 
matter and quintessence, and ignore the neutrinos as 
well as the distinction between baryons and CDM. Each 
component i (i = r, to, Q) has an energy density pi and 
pressure pi . The sum of the energy densities determines 
the Hubble parameter via the usual Einstein equation, 
3H 2 = p = J2iPi- The equation of state = Pi/pi is 
1/3 and zero in the case of radiation and matter, respec- 
tively. The background energy density and the pressure 
of the quintessence field are: 



pq = \q 2 + v{Q) 



PQ = \Q 2 - V{Q) , (2) 



where the background scalar field Q(t) obeys the equa- 
tion Q + 3HQ + Vq = 0. This means that wq is in 
general time-dependent. 



a 
cT 



i ■ 

0.75- 
0.5- 
0.25- 
- 
-0 .25 
-0.5- 
-0.75 
-1 - 



1 

. 75 
0.5 
.25 



-0 . 25 
-0.5 
-0 . 75 

-1 



10 



10" 



10 



1 



1+z 



FIG. 1. Densities of radiation (solid line, black) and matter 
(dotted line, blue), and the equation of state for the scalar 
field (dashed line, green), as a function of redshift for a model 
with potential V(Q) = M i e f/Q . In this plot M 4 = 10 
and / = 1 M p i (with 8ttG = A-/" 2 .) 



The main requirement of quintessence is that it starts 
to dominate the energy density of the Universe only at 
the present time, with an equation of state wq(z = 0) = 
wqq < 0. Conservative phenomenology dictates that 
PQo/po = ^qo < 0.8 — to allow time for galaxy for- 
mation — and that wqq < —0.5 — to accommodate the 
SNIa data j2(|. In addition, nucleosynthesis demands 
that fl Q < 0.2 at z ~ 10 9 . 

Model building should obtain these values and still 
manage to solve (or at least alleviate) the coincidence 
problem without too much fine tuning (see, e.g., p7|.) 
The main problem seems to be that in order to solve 
the coincidence problem one needs a period of tracking, 
but it is hard to get a period of tracking and still ob- 
tain wqq < —0.5. It is certainly possible to construct 
potentials which implement both conditions |7j , however 
we prefer to look at simpler potentials which contain fea- 
tures which are generic to most of the models. As it turns 
out, the relevant features of the cosmological perturba- 
tions do not depend on the specific form of the potential, 
but only on its generic phenomenology. 

In a typical scenario, nicely reviewed in the scalar 
field starts out subdominant deep in the radiation era, 
in a kinetic phase with an equation of state wq = +1 
(see Fig. 1.) The kinetic energy i(j 2 eventually decays, 
leaving only the quintessence potential energy, which is 
nearly constant and as a result wq — > —1 (of course, 
the quintessence field could also start already in the 
potential-energy dominated regime) . This is the so called 
potential phase. When Qq becomes of order unity, the 
quintessence field undergoes a transition which puts it 
into a tracking regime, where it follows approximately 
the equation of state of the background. Finally, at some 
point late in the matter era Q starts to dominate and the 
Universe begins the accelerated expansion phase that we 
observe today (the Q-dominated phase.) 

As long as we keep away from the time of equal matter 
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and radiation (z eq ~ 10 4 in our flat models with Hq = 65 
Km s -1 Mpc -1 ), one of the two barotropic fluids (radia- 
tion or matter) can be neglected. The equation of state 
of the total matter content then reads: 



W = - = Wf + (w Q - w f )Qq 



(3) 



where the subscript F stands for either radiation, when 
t <C t eq , or matter, when t 3> t eq . The total speed of 
sound c 2 has a simple expression as well: 



2 _ P Q , 2 \ 
c s = ~ = w F -\ — (c Q - w F ) , 

P P + P 

where Cq = Pq/ Pq- It is useful to note that 
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When Cq ^constant, then Q oc Vq/H. Taking the time 
derivative of c?, we also obtain: 
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When wq is approximately constant, then from Eq. (|^) 
Cq is constant as well. By using the background relation 
H/H 2 ~ -3(1 + w F )/2 we obtain 



Vqq ~ -H 2 (l - c 2 Q )(w F +c 2 q + 2) = aH 2 



(7) 



We observe that quintessence fluctuations are effectively 
massless {Vqq ~ 0) during the kinetic phase (cq = 1) 
and the potential phase (cq = — 2 — w F 0). 

A good quantity to measure how closely the 
quintessence field Q tracks the background is the quantity 
7 = VqqV/(Vq) — which was named r in ||. When 
wq is approximately constant in time, then: 
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Wp - Wq 

2{1 + wq) 



(8) 



If tracking is exact (as is the case in the exponential po- 
tential models), w F = wq, then 7=1 and Vqq/Vq — 
Vq/V. _ 

If wq is approximately constant, then the kinetic and 
potential energies of the scalar field Q must be pro- 
portional to each other. Taking the time derivative of 
Pq oc Q 2 oc V we obtain that during a tracking phase: 

V q , H 

= —3(1 + wq)— ~ constant , (9) 

where we have used the fact that Q 2 oc a~ 3 ( 1+WQ ^ to find 
H/Q oc a - 3(wF - w ^ oc constant, if w Q ~ w F . 



We note that during the phases in which wq is constant 
Hq goes as: 



= — a 3(ti>TOT-tUQ) 

Ptot 



(10) 



Therefore in the radiation epoch VLq redshifts as a~ 2 
during its kinetic phase, and it grows as a 4 during its 
potential phase. During the tracking regime Hq is ap- 
proximately constant (depending on how accurate is the 
tracking) and in the Q-dominated regime it is evidently 
constant (since £Iq — > 0(1) when Q dominates.) 



III. EVOLUTION OF PERTURBATIONS 

We compute the evolution of the cosmological pertur- 
bations in longitudinal (or conformal-Newtonian) gauge 
PI: 



ds 2 = (1 + 2$)dt 2 - a 2 (t)(l - 2$)dx : 



(11) 



In the matter sector, the perturbed energy density and 
pressure of quintessence are, as usual, 



Spq = Q5Q - Q 2 $ + VqSQ 
= QSQ - Q 2 $ - VqSQ 



(12) 
(13) 



where the Fourier transforms of the scalar field fluctua- 
tions obeys the equation: 

... A; 2 
SQ + 3H5Q + -^SQ + VqqSQ = +4Q$ - 2V Q <S> . (14) 

The density fluctuations in radiation and matter, on the 
other hand, obey the conservation equations for the den- 
sity contrasts S F = Sp F /p F : 

i F -3{l + w F )^ = (l + w F )-^rV F , (15) 

a{t) 

where V F is the fluid velocity. In the long wavelength 
limit (k <C aH) this equation is extremely useful, since 
it reads: 

S F - 3(1 + w F )$ = constant = 8 F - 3(1 + w F )^ , (16) 

where the superscript i indicates that the fluctuations 
have been evaluated at some initial time tj. We can also 
combine Eq. ( |l5| ) for 5 F with the equation for conserva- 
tion of momentum, 



V F =*(-^6 F + *) , 
a(t) \1 + w F ) 

and obtain a simple second order equation for S F : 

k 2 N 



(17) 



5 F " + w F k 2 6 F = 3(1 + w f ) [ $ " - y $ J , (18) 
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where a prime denotes, as usual, derivative with respect 
to conformal time 77, where d/dr) = a(t)d/dt. 

We consider now the evolution of the quintessence field 
perturbations. At first, let us neglect the metric pertur- 
bations, i.e. take Eq. ( fl4| ) without its right-hand side. 
By using the rescaled variables SQ = a x ^ 2 SQ, Eq. (|l4| ) 
can be rewritten as: 



SQ" + H8Q' 



k 2 + a 2 Vno~^- — 
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2a 
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SQ = 



(19) 



where TL = a 1 /a = aH. 

In a radiation dominated universe a oc 77 and the term 
proportional to the second time derivative of the scale 
factor vanishes; when Eq. (0) holds, then the solutions 
for the field perturbations in rigid space-times are: 



5Q~n 



-1/2 



where 



--a. 



(20) 



If a ^ both the solutions decay in time. If a — > then 
there is a constant mode. 

The argument above holds for a matter dominated uni- 
verse as well: if a = there is a constant mode, oth- 
erwise both the solutions decay in time (for a matter 
dominated universe the appropriate rescaled variable is 
SQ = a 3 / 4 SQ.) 

The inclusion of gravitational fluctuations in Eq. fll4| ) 
leads to a constant solution for the quintessence pertur- 
bations in the long- wavelength limit. From Eq. (|l4|) we 
immediately see that in this limit there are constant so- 
lutions $(t) $ c and SQ(t) -> 5Q C : 



SQ C 



-2^<& c 



(21) 



as long as Vq/V^qq is approximately constant. But this 
is precisely what happens during the tracking regime: 
from Eqs. (0) and (0) we see that Vq/Vqq ~ V/Vq ~ 
constant in the tracking period. We stress that the type 
of solution (^l|) does not hold in the kinetic and potential 
phases, since in these cases Vqq — ► 0. 

We can use the — component of the Einstein equa- 
tions to relate the Newtonian potential to the energy den- 
sities of other fluids 



6H 2 $ - 6H$ - 2—$ = Sp r + 5 Pm + S P Q 



(22) 



In the long wavelength limit, assuming that $ is station- 
ary and ignoring the subdominant barotropic fluid we 
obtain: 



2$ ~ n F s 



F°F 



flnS, 



Q°Q 



(23) 



In the regime described by the attractor in Eq. ( |2l| ) 
the perturbed energy density for quintessence, defined in 
Eq. ([l2]), reduces to: 
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(24) 



and the perturbed pressure defined in Eq. (O) is: 



-2p Q $ + 2$(Q 2 -!/) 



7-1 

7 



+ 0( 7 ). (25) 



Using now the background identities 3H 2 = p and ~ 
1 — Hp into Eq. (|24|), we obtain the density contrasts 
as functions of the Newtonian potential in the tracking 
regime: 



Q 



-2$ c 



(26) 



Notice that during tracking, quintessence and the dom- 
inant fluid species are in effect indistinguishable (Sq ~ 
S F ), consequently we expect isocurvature perturbations 
to be suppressed during that period. 

A solution corresponding to (21) — though in the syn- 
chronous gauge — was first obtained in the case of expo- 
nential potentials, for which tracking is exact An 
attractor for quintessence perturbations has also been 
conjectured in 0. As we have shown, the approximate 
solutions (|l|) and ( p6| ) hold for any potential, just as 
long as there is tracking, and, of course, they exist in any 
gauge. 

We can also combine Eqs. ( p6| ) and ( [l6| ) to obtain a 
relationship between the initial and final values of the 
Newtonian potential: 



-- (5 i 
6 1 r 



4<T 



(27) 



where 5 l r and <I> 1 are initial conditions for the radiation 
density contrast and the Newtonian potential respec- 
tively. Notice that these initial conditions can be speci- 
fied even at a time when quintessence is dominating, and 
even if the radiation contrast and the Newtonian poten- 
tial are initially not constant. 

We have numerically verified formulas (|2l]), (|24|), ( |2"6| ) 
and (|27]) for several scenarios and initial conditions. Take 
for example the scenario whose background appeared in 
our Fig. 1. Two typical initial conditions for the field 
perturbations are plotted in Figs. 2A and 2B (solid 
and dashed lines), together with the attractor solutions 
(thin red lines): as the scalar field starts to approach 
the tracking regime at z ~ 10 12 , the perturbations start 
to converge around the attractor solution. As seen in 
Fig. 2A for the field perturbations SQ, the attractor of 
Eq. ( ^l| ) is a very good approximation even after the 
tracking phase has ended — that is, Eq. (pll) is a very 
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FIG. 2. Quintessence field fluctuations SQ (panel A) and 
quintessence energy density contrasts Sq (panel B) for two 
different initial conditions for SQ and SQ (solid and dashed 
lines.) We have chosen initial conditions such that <E> ~ 10~ 5 
at Zd C c = 10 3 . In panel A the quintessence fluctuations 
(dashed and solid black lines) are compared with the attrac- 
tor — 2§V,q /V,qq (thin red line). Notice that even after the 
tracking phase ends the analytic approximation of Eq. ( |2~i] ) re- 
mains very good. In panel B the quintessence energy density 
contrasts (solid and dashed black lines) are compared with 
— 2<E> (thin red line), verifying the second approximation, Eq. 
(^6j). The approximation becomes worse as the tracking ends. 
In both plots the wavelengths of the modes cross the Hubble 
radius at z = 0. 



good approximation during the quintessence-domination 
period as well, even though SQ is not constant anymore. 
Fig. 2B shows how the quintessence density contrasts 
converge to —2$. In fact, for a wide range of initial 
values the scalar field perturbations end up at the same 
solution 5Q C after tracking. 

This can be also seen in Fig. 3, which is the phase di- 
agram for the perturbations with different initial condi- 
tions shown in Fig. 2. During tracking the solutions spi- 
ral down to the attractor (solid and dashed lines.) When 
tracking ends the attractor disappears, but by that time 
most modes have settled down to the same value, and 
their evolution is henceforth the same (see the thin red 
line in Fig. 3 which springs from the attractor point.) 
We note that this attractor occurs generally only during 
tracking and for long wavelengths: when the wavelength 
becomes important (of the order of the inverse effective 




FIG. 3. Phase diagram for quintessence fluctuations 
with two different initial conditions. The attractor point 
( — 2&V t Q/V t QQ, 0) is a transient attractor, valid only for long 
wavelength modes. As soon as the tracking regime of the 
field Q terminates, the attractor disappears and the field per- 
turbations start to evolve (green straight line.) However, by 
the time tracking is over, most long-wavelength solutions have 
already converged to the same value, and thereafter their evo- 
lution is almost indistinguishable. 

mass []), there is a sensitivity with respect to the initial 
conditions of the quintessence perturbations. 

Summarizing the results of this section: in the kinetic 
and potential phases there is a constant long wavelength 
solution which is a linear combination of the initial field 
fluctuations and of the gravitational potential, as one can 
see from the initial evolution in Fig. 2A. In the track- 
ing period the quintessence perturbation stabilizes at the 



attractor solution (21) and it remains at that constant 
value until the perturbation reenters the Hubble radius, 
or until quintessence starts to dominate the background. 



IV. EVOLUTION OF ISOCURVATURE 
PERTURBATIONS 

Cosmological fluctuations are often characterized in 
terms of the gauge-invariant curvature perturbation on 
comoving hypersurfaces, defined as pfl: 



2$/H + <$> 
^ = 3 1 + w 



(28) 



The time variation of the intrinsic curvature £ is given, 
on large scales, by the non-adiabatic pressure or, equiva- 
lently, by the amplitude of the isocurvature perturbations 
— see Eq. ([!]) . Therefore, if the non-adiabatic pressure 
vanishes then ( is constant. 



x If the effective mass is of the order of the Hubble radius, 
this occurs when the perturbations reenter inside the Hub- 
ble radius. It occurs at the same time if the effective mass 
vanishes (as happens during the kinetic and potential phases). 
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The pressure perturbations can be split into adiabatic 
and non-adiabatic components: 



Sp 



Sp 



5p + 



5r=o 



Sp 

sr 



ST, 



(29) 



<5p=0 



where we define the entropy perturbation <5r as 28 
dp Sp 



ST 



P 



P 



(30) 



As a consequence we have that the adiabatic pressure 
perturbation is given by 



<5pad = c 2 s Sp : 



Sp 



Sp = ^Sp , 



(31) 



sr=o 



and the non-adiabatic pressure is given by the second 
term in the right-hand side of the definition (j29|): 



Sp 



nad 



Sp 

ST 



ST = pST 



(32) 



s P =o 
Sp - c 2 s Sp 



E 



(Spi - c 2 s Spi 



We now give an analytic description of the time evolu- 
tion of isocurvature perturbations for long wavelengths. 
We work perturbatively assuming £ constant, and com- 
pute the non adiabatic pressure pST: when the inte- 
grated effect is large, it means that isocurvature pertur- 
bations cannot be neglected. The results of this analysis 
confirm the methods of the previous section (where $ 
was assumed constant) and are in agreement with the 
numerical analysis. 

Using Eq. (J3^) we find that the non-adiabatic pressure 
is given by: 



<5Pnad = {w r - C 2 S ) Sp r + Sp Q 



By using Eqs. ( 
that in general 

Sp n 



ad 



p + p 



and (0) into Eqs. 



= o(n Q ) x (s r + s Q ) . 



(33) 
one can check 

(34) 



Therefore, when the quintessence contribution to the to- 
tal energy density is very subdominant, the isocurvature 
contribution is small. However, since the isocurvature 
contribution to £ is an integrated effect, £ ~ tSp n &d, we 
should study the time evolution of each term which enters 
into the definition of the non-adiabatic pressure (|33J) . 

The first term in the right-hand side of Eq. (|33|) is 
proportional to p r S r and leads at most to a logarithmic 
increase of (. 

Among the contributions from the pressure and the 
energy density of the scalar field we neglect the first term 
in the right-hand sides of both Eqs. (12) and (|l3|). The 
term Q 2 <& is suppressed during the kinetic phase, but 
leads to a growth as a 7+3l " TOT in the left-hand-side of Eq. 



(p4|) in the potential regime. During the tracking regime 
Q 2 gives approximately a constant contribution to (|34|). 
The term V,q SQ depends explicitly on the quintessence 
fluctuations: it decays during the kinetic phase and it 
grows less rapidly than Q 2 $ during the potential regime. 
However, it leads to a growth as a 3 / 2 ( 1 +™TOT) d urm g the 
tracking regime. 

It is clear from Eq. (^) that for two barotropic flu- 
ids with the same equation of state and the same den- 
sity contrast the non-adiabatic pressure should be zero. 
Therefore for exact tracking, quintessence and radiation 
equilibrate to give zero non-adiabatic pressure. We can 
also compute the non-adiabatic pressure in the tracking 
regime by using Eqs. (|33|), (4) and (||): 



Sp 



nad 



Q 2 



(w r - Cq) (p r S r - 2$p Q ) 



P+P 
+ 0[( 1 -l)p Q ^ 



(35) 



As expected, the non-adiabatic pressure vanishes for ex- 
act tracking (7 = 1), otherwise it is small, but not van- 
ishing, during tracking. In general <5p na d is proportional 
to /9q, therefore small in many models. 

There is however still one possibility that allows for 
significant isocurvature fluctuations from quintessence: 
this happens when the tracking phase starts only at a 
relatively late redshift, z ~ 10 5 — 10 3 (see also the next 
Section, Figs. 4, 5 and 6.) During the transition from 
the potential phase to the tracking phase, we can have 
both £Iq non-negligible, and the equation of state and 
speed of sound of quintessence differ substantially from 
those of the background fluid. If the quintessence den- 
sity contrast is of the same order of the matter density 
contrast at decoupling time z^cc ~ 10 3 , the isocurvature 
perturbations can leave an imprint on the CMBR. Later, 
as tracking forces the field perturbations to the attractor, 
the isocurvature fluctuations are temporarily depressed, 
at least until Q starts to dominate at z ~ 0. 



V. INITIAL CONDITIONS: ADIABATIC OR 
MIXED? 

In this section we address the issue of how initial con- 
ditions of quintessence perturbations can be used as tools 
in CMBFAST Q . These initial conditions are set after 
nucleosynthesis, at z ~ 10 9 . In most of the literature, the 
initial conditions for the quintessence fluctuations are set 
up by requiring adiabaticity with the other components 
However, the notion of a purely adiabatic pertur- 
bation (as well for a purely isocurvature one) is an in- 
stantaneous notion for a multifluid system. Moreover, 
because of the unthermalized nature of quintessence, the 
adiabatic condition for this component is even less justi- 
fied. 

The adiabatic condition ||D is usually defined as the 
vanishing of the relative entropy and its time derivative: 



G 



SrQ = , 
SrQ = , 

which reduces, in longitudinal gauge, to: 

-l 



SQ = Q 2 
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The relative entropy between the radiation and the 
quintessence components is defined asQ 



S, 



tQ 



S r 



1 



1+WQ 



(40) 



From the above relation we immediately understand that 
S r Q — in the case of exact tracking. Indeed, two fluids 
with the same equation of state and the same density 
contrast are undistinguishable. For long wavelengths re- 
lations (p^)-(B9) reduce to: 



SQ ~ -^-S r * — $ 
^ AH 2H 



SQ~Q 



S r 

T 



HQ J 



-§•<> 



HQ 



(41) 



(42) 



where the second line of Eq. (|42j) holds only if radiation 
is the dominant component. 

The conditions ©-© should be compared with 
the solutions during the kinetic and potential phases, 
SQ =constant and SQ = 0, or with the attractor solu- 
tion Eq. (pTi) which applies during the tracking phase: 




VnV. 



QQQ 



V 2 
V ,QQ 



(43) 
(44) 



It is not hard to see from Eq. (|4) that in the track- 
ing regime SQ -> 0, since V Q /V QQ ~ Vqq/Vqqq. 



2 Notice that the definitio n (MOI) of entropy is different from 
the SF introduced in Eq. (BOh^We consider this more stan- 
dard definition S r Q as well since it is this quantity which is 
used in much of the literature to define the adiabatic condi- 
tions, and we want to show directly the difference between 
the initial conditions (^)-(^) and the ones defined by the 
attractor or by some previous dynamics. 




1+z 

FIG. 4. Background evolution in a quintessence scenario 
with the Ratra-Peebles potential V(Q) = M 4+a Q~ a 1 where 
a — 6 and M — 10 -6 ' 15 . Plotted are the densities of radiation 
(solid line, black) and matter (dotted line, blue), and the 
equation of state for the scalar field (dashed line, green), as a 
function of redshift. 



Using this fact together with Eq. (||), we also obtain 
that Eqs. ( pi] ) and (|i^) are similar during the track- 
ing regime. Therefore, as expected, if the scalar field is 
tracking normal matter, then the adiabatic condition is 
approximately satisfied. 

We illustrate the previous discussion in Figs. 4, 5 and 
6. The background model is plotted in Fig. 4, and in 
Figs. 5-6 we present possible scenarios for the pertur- 
bations, as well as a comparison with the cosmological 
perturbations in the ACDM case. 

The upper (thin red) curves of Fig. 5 plot the gauge- 
invariant curvature £, defined in Eq. (|28|). The lower 
(black) curves are the Newtonian potential <I>. All plots 
in Fig. 5 have been normalized so that $ = 5.10~ 6 at z = 
10 9 , and the wavelengths corresponds to modes which are 
crossing the Hubble radius at the present time (z = 0.) 

The solid lines of Fig. 5 correspond to the cosmological 
perturbations of a fiducial ACDM scenario with adiabatic 
initial conditions. Notice that, as usual, since £ remains 
constant, <& has to change by 9/10 after z eq — 10 4 . 

The long-dashed lines of Fig. 5 are the perturbations 
in a scenario (ICDM) where, in addition to the adiabatic 
mode, CDM and radiation have an initial relative isocur- 
vature of S r CDM — 3<I>/2. 

The short-dashed lines are the perturbations in the 
case of adiabatic initial conditions (AIC) between all 
components at z ~ 10 9 . The dotted lines are the per- 
turbations in the case where we choose zero isocurvaturc 
between radiation and CDM, and SQ = 10~ 3 , SQ = 
initially (QIC.) This last set of initial conditions (QIC) 
is motivated by the fact that the field perturbations are 
constant during the kinetic and the potential phases — 
see, e.g., Fig. 2A. 

Notice the identical late isocurvature effect in AIC and 
QIC. The signal of this effect is the extra growth of the 
Bardeen parameter £ at late times (compare the dotted 
and short-dashed lines with the solid line in the interval 
< z < 10 in Fig. 5.) This effect is independent from 
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1+z 

FIG. 5. Newtonian potential (lower, black curves) and 
gauge- invariant curvature £ (upper, red curves), normalized 
by the condition that $ = 5.1CT 6 at z = 10 9 . The solid lines 
are a fiducial ACDM model. The remaining lines correspond 
to different initial conditions for the fluctuations of the radia- 
tion, matter and quintessence components. The long-dashed 
lines correspond to the case where there is no isocurvature 
component between Q and radiation, but there is an initial 
isocurvature component between radiation and CDM. The 
short-dashed line corresponds to the case of pure adiabatic 
initial conditions (AIC). The dotted lines correspond to the 
case of SQ = 10 -3 , SQ = initially (QIC). All modes cross 
the Hubble radius at z = 0. 



- Vio 

-6 -4 -2 1 

10 10 10 

1 + Z 

FIG. 6. Newtonian potential (lower, black curves) and 
gauge-invariant curvature £ (upper, red curves), for a wave- 
length which crosses the Hubble radius at decoupling time 
between matter and radiation. The normalizations and the 
initial conditions are the same of Fig. 5. As in Fig. 5, the 
short-dashed lines correspond to the case of pure adiabatic 
initial conditions and the dotted lines correspond to the case 
of SQ = 10~ 3 , SQ = initially. Even for smaller scales, the 
effect of an isocurvature component is a transient effect. For 
wavelengths larger than the Hubble radius, the attractor was 
responsible for explaining the transiency, while for smaller 
scales the explanation is the decay of quintessence fluctua- 
tions inside the Hubble radius. 



the initial conditions for the quintessence fluctuations. 
In fact, as already emphasized, the notion of adiabaticity 
(as well as pure isocurvature) is an instantaneous one, 
imposed at an initial time, and it does not persist in a 
multi-fluid system. We stress that this effect is distinct 
from the change of the Newtonian potential $ which is 
due to the late change in the equation of state of the back- 
ground, which can be seen in pure form in the ACDM 
case (solid line in Fig. 5.) A similar change in 4> occurs 
also at z oq — 10 4 , while £ remains constant across the 
transition between radiation- and matter-domination. 

Notice also the early isocurvature effect in the QIC 
case (oscillations of the dotted lines in Fig. 5.) This 
means that there is a substantial non-adiabatic pressure, 
and hence a large isocurvature perturbation, in this sce- 
nario. As we discussed at the end of the previous Section, 
the reason why isocurvature fluctuations can become im- 
portant is that the "tracking regime" of the background 
model only starts relatively late, at z ~ 10 5 . But the 
quintessence field perturbations need some time to con- 
verge to the attractor solution. During this time £Iq 
becomes non-negligible, and since the Q component still 
behaves quite differently from the dominant background 
fluid, there can be substantial isocurvature perturbations 
for wavelengths which cross the Hubble radius at z ~ 0. 
A similar transient effect occurs also for smaller scales, 
as shown in Fig. 6, for a wavelength which crosses the 
Hubble radius at the decoupling time. An isocurvature 
component from the quintessence field leads only to a 
transient effect, since quintessence fluctuations converge 
to an attractor for long wavelengths, and decay in time 



inside the Hubble radius. However, this early transient 
isocurvature effect can leave an imprint in the spectrum 
of the CMBR anisotropies. 

VI. DISCUSSION AND CONCLUSIONS 

We have studied the evolution of cosmological pertur- 
bations in quintessence models, with particular attention 
paid to isocurvature modes generated by quintessence. 
We have shown that these isocurvature modes are 
generic, with the exception of the case in which 
quintessence mimics exactly radiation. This occurs in 
the tracking phase of models with exponential potentials 
However, these models are unappealing from the 
phenomenological point of view because their equation 
of state for quintessence is not negative at the present 
time. 

When tracking is not exact, then isocurvature modes 
are non-vanishing. We found an attractor solution for 
long wavelength quintessence fluctuations in the track- 
ing regime. This allows an estimation of the amount 
of isocurvature fluctuations in the tracking regime for 
any quintessence model that display a tracking period. 
In the other phases which usually occur in models with 
quintessence — the kinetic and potential phases — the 
quintessence fluctuations have also a constant mode, 
whose actual value is determined by the previous evo- 
lution (including inflation.) The contribution of isocur- 
vature fluctuation to the adiabatic mode grows during 
the potential, tracking and Q-domination phases. 



We have discussed the assumption of adiabaticity in 
the context of the setting of initial conditions used in 
numerical codes such as CMBFAST. As already empha- 
sized, because of the lack of thermal equilibrium, this 
assumption is not justified for quintessence. However, in 
the tracking case, the conditions set by the attractor so- 
lution are in most cases quantitatively indistinguishable 
from the adiabatic conditions. Indeed, tracking seems a 
gravitational mechanism — alternative to thermal equi- 
librium — which tends to reduce isocurvature modes be- 
tween quintessence and the background fluid. For mod- 
els with a tracking phase which starts early in time, we 
expect a weak dependence on the initial conditions for 
the quintessence fluctuations. In the stages prior to the 
tracking phase the "initial conditions" (understood as the 
values of the held perturbation and its time derivative at 
a redshift of z ~ 10 9 ) depend on the conditions set by 
inflation, and could be different from the adiabatic ones. 

We have identified a late isocurvature effect for long 
wavelengths due to quintessence. Of course, this is due 
to the fact that Q dominates at late times, and it is qual- 
itatively independent from the model considered. There- 
fore, the theoretical explanation of the long wavelength 
evolution of the Newtonian potential in Q models is a su- 
perposition of two effects: the change of the equation of 
state and the growth of Q on super-Hubble scales. This 
late isocurvature effect could be useful in order to dis- 
tinguish a cosmological constant model from the models 
with a quintessence component. 

The observational relevance of isocurvature modes gen- 
erated by quintessence is weakened by the decay in time 
of quintessence perturbations inside the Hubble radius. 
For this reason the isocurvature mode in the Q-radiation 
sector are very different from those in CDM-radiation. 
This effect was also very appealing in order to minimize 
the effects of the inclusion of this extra component on 
structure formation. In the models examined here, the 
suppression of Qq during the kinetic phase play also a 
crucial role in order to weaken the effect of some initial 
isocurvature modes generated by quintessence. Even if 
the upper bound during nucleosynthesis for Q,q is ~ 0.2 , 
the kinetic regime suppresses FIq down to 10~ 15 — 10~ 20 
in the models which we have analyzed. Therefore our 
considerations may be more relevant for models in which 
Qq is closer to the upper limit, as, for instance, in the 
models with a modified exponential potential j3fj] . 

It is therefore interesting to study the impact of isocur- 
vature fluctuations in quintessence models using numer- 
ical tools such as CMBFAST ]3(| in which quintessence 
perturbations are included. In particular, it is possible 
to construct models where FIq is not so suppressed or 
in which there is a late tracking phase, in which case 
isocurvature modes generated by quintessence could lead 
to observable effects in the temperature anisotropies of 
the cosmic background radiation. 
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VII. APPENDIX A 

Here we present the equations that were evolved nu- 
merically. Combining the — and i — i Einstein field 
equations we can eliminate the radiation perturbations, 
and obtain the following equation for the quintessence 
and matter perturbations (as explained in the text, we 
work with units such that 8irG = 1, and we have ignored 
the distinction between baryons and CDM): 

3$ + 15H<f> + (l2H 2 + ^) $ (45) 

= -\s mPm + QSQ - Q 2 $ - 2V Q SQ . 

The equation of motion for the scalar field perturbations 
is given in Eq. (|l4|). The matter density contrast obeys 
the energy conservation equation: 

S m - 3$ - -v m = , (46) 
a 

where v m is the matter fluid velocity, which itself satisfies 
the momentum conservation equation: 

k 

Vra +Hv m + = . (47) 

a 



VIII. APPENDIX B 

Here we present the evolution of pure isocurvature 
modes generated by the quintessence sector. Consider 
contributions from quintessence and another component 
X to the metric perturbation, which cause a vanishing 
contribution to $ and 4> |l0| at some fixed initial time 
in the energy and momentum constraint. This requires 
that: 

n Q s Q = -n x s x (48) 

QSQ = --v x px, (49) 

which identifies an isocurvature density mode ( f48| ) and 
an isocurvature velocity mode (^) |19[] . The third com- 
ponent Y is the only one which leads to a curvature per- 
turbation (for long wavelengths —2$ = TySy, as follows 



9 




10 

1 

1+z 



FIG. 7. Gauge-invariant curvature £ for Q-CDM isocurva- 
ture (solid lines) and Q-radiation (dashed lines) initial con- 
ditions for a wavelength which crosses the Hubble radius at 
z — 0. The normalization is made in such a way that the 
amplitude of the Newtonian potential is roughly 10 -5 when 
the growth saturates. 



from Eq. (|23|).) The components X and Y are related 
by adiabaticity (Sxy — 0.) Therefore, the quintessence 
density contrast at fixed initial time is: 



ox 



n 



X 



2$ 



1 + Wx 

1 



Wy 



X 



(50) 



In Fig. 7 the evolution of the Bardeen parameter £ 
for initial pure Q-CDM (X = CDM, solid line) and Q- 
radiation (X = radiation, dashed line) isocurvature ini- 
tial conditions at z = 10 9 for the model of Fig. 4. At the 
initial time the Newtonian potential is so small in order 
to ensure 5q < 1. However, in both cases - but in par- 
ticular for the Q-radiation isocurvature initial conditions 
-, Eqs. force Q fluctuations to very large val- 

ues, either initially and in the subsequent evolution. The 
evolution of Q fluctuations for these type of initial con- 
ditions is not described by the attractor (|Tj) - i.e. the 
inhomogeneous solution to the differential equation for 
the quintessence fluctuations -, but by the homogeneous 
solution to Eq. (|l4j). Indeed, in Fig. 7 the growth of 
the Bardeen parameter saturates roughly when tracking 
starts in Fig. 4, i.e. when the homogeneous solution for 
quintessence fluctuations start to decay, as described by 
Eqs. (19)-(20) with a ^ 0. 
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